Polarization properties in the transition from below to above lasing threshold in 
broad-area vertical-cavity surface-emitting lasers 
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For highly divergent emission of broad-area vertical- cavity surface-emitting lasers (VCSELs) a ro- 
tation of the polarization direction by up to 90 degrees occurs when the pump rate approaches the 
lasing threshold. Well below threshold the polarization is parallel to the direction of the transverse 
wave vector and is determined by the transmissive properties of the Bragg reflectors that form the 
cavity mirrors. In contrast, near-threshold and above-threshold emission is more affected by the 
reflective properties of the reflectors and is predominantly perpendicular to the direction of trans- 
verse wave vectors. Two qualitatively different types of polarization transition are demonstrated: an 
abrupt transition, where the light polarization vanishes at the point of the transition, and a smooth 
one, where it is significantly nonzero during the transition. 

PACS numbers: 42.55.Px, 42.60.Jf, 42.25.Ja, 05.40.-a 



Oh 

d 



> 

(N 

oo 
m 

On 
O 

a\ 

O 



I. INTRODUCTION 

In the last decades vertical cavity surface emitting 
lasers (VCSELs) have played an increasing role in sci- 
entific research and applications [1]. One of the features 
of VCSEL design is the possibility to obtain large two- 
dimensional apertures which are quite homogeneous and 
have a small polarization anisotropy. 

The polarization behavior in such lasers resulting from 
the competition of stimulated and spontaneous emission 
has been a subject of many investigations [H, 0, 0, [H, 
It is known that the polarization properties of small and 
medium aperture VCSELs above flii[iO,[iil,[il and 
below 0, M\ threshold are determined mainly by the in- 
tracavity anisotropics. Well below threshold the polar- 
ization degree is reduced dramatically, in VCSELs [4] as 
well as in edge-emitting semiconductor lasers [13.] (which 
possess much higher intracavity anisotropics). However, 
below threshold the polarization coincides with the at- 
threshold one. 

With increasing size of the aperture, off-axis emission 
becomes important. Because the cavity resonance is dif- 
ferent for different transverse modes, the ones with the 
best alignment between the cavity resonance and the 
gain maximum of the active medium have the largest 
gain [ii, [iH. Recently it was shown [11, [13, [3 that 
for strongly off-axis emission the intra-cavity anisotropics 
play only an auxiliary role in polarization selection above 
threshold. In contrast, the polarization-selective prop- 
erties of reflection and transmission of the distributed 
Bragg reflectors (DBRs) forming the cavity mirrors are 
much more important. Above threshold, the DBR TE 
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modes (which are in paraxial approximation perpendic- 
ular to the transverse component of the wave vector, i.e. 
s-waves) have higher reflectivity^ and the cavity quality 
factor for the TE modes is larger than for the TM modes. 
This determines the polarization of the above threshold 
emission, which has an overall tendency to be perpendic- 
ular to the transverse wave vector ( "90-degree rule" ) [17| . 
Later investigations established that the above-threshold 
polarization is also strongly affected by the transverse 
cavity boundaries, i.e. the waveguide formed by the oxide 
confinement [18]. This leads to strong deviations of the 
above-threshold polarization state from the "90-degree 
rule" for transverse wave vectors with directions not par- 
allel to either of the device boundaries [18]. The data in 
[l^ (for lasers with a square aperture) and in [l^ (for 
lasers with a circular aperture) indicate also that the po- 
larization direction for off-axis light is different below and 
above threshold but there is no detailed investigation. 

In this work we characterize the polarization properties 
of off-axis below-threshold emission and show that the 
polarization direction is governed mainly by the trans- 
missive properties of the top DBR. The transmissivity is 
larger for the TM Bragg modes (parallel to the transverse 
wave vector, p- waves) than for the s-waves, resulting in a 
"0-degree rule" for polarization selection, i.e. the polar- 
ization is parallel to the transverse wave vector. 

We consider, both theoretically and experimentally, 
the transition from the below-threshold to the above- 
threshold polarization state for highly divergent VCSEL 
emission. The nature of the transition depends critically 
on the orientation of the polarization of the final (lasing) 
state. When the final polarization obeys the "90-degree 
rule" and is perpendicular to the polarization state well 
below threshold, the transition is very abrupt and the 
light is unpolarized at the point of transition. On the 
other hand, when the polarization of final state is not 
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orthogonal to the initial one, the transition is consider- 
ably smoother, and the emission retains a relatively large 
degree of polarization during the transition. Our theory 
predicts also that far below threshold the main principal 
axes of the intra-cavity and extra-cavity field are per- 
pendicular to each other due to the strong anisotropic 
filtering of the light coupled out via the DBRs. 

In the next section we describe the experimental setup; 
in Sect. IIIII the experimental results are reported. In 
Sect. HVl a theoretical model for the description of the 
transition is developed, analyzed and the results com- 
pared to the experimental observations. Concluding re- 
marks are in Sect. \V\ 
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EXPERIMENTAL SETUP AND METHODS 
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FIG. 1: (Color online) Experimental setup. The VCSEL is 
set into an air-tight box to avoid condensation water. Polar- 
ization optics: half-wave plate and linear polarizer. 



The VCSELs under study are oxide-confined top- 
emitters with a square aperture of 40 x 40 jitm^ that 
are packaged in TO-type housings without caps. The 
emission wavelength is around 780 nm. The lasers con- 
sist of two highly reflective DBRs (top mirror: 31 lay- 
ers, bottom mirror: 47 layers) with three 8 nm thick 
Alo.11Gao.89As quantum wells in between. Together with 
several Alo.36Gao.64As spacer layers and the GaAs sub- 
strate the whole structure is about 10 \xm long. In order 
to reduce the electrical resistance of the lasers the inter- 
faces between different semiconductor layers are graded. 
A laterally oxidized layer above the active region provides 
current and optical confinement. 

The devices are electrically pumped with a low-noise 
DC current source between and 30 mA. The typical 
lasing threshold at °C heat sink temperature is 15 mA. 
The VCSEL is mounted on a copper plate that is at- 
tached to a thermo-electric cooling element enabling tem- 
perature control between 40 °C and —35 °C. We remark 
that the actual device temperature is strongly influenced 
by the driving current due to Joule heating effects. Tem- 
perature values given here are for the heat sink only. The 
device temperature can be inferred indirectly from opti- 
cal spectra and changes with current by a factor of about 
0.9 K/mA. The device temperature changes the detun- 
ing between the cavity resonance and the maximum gain 
frequency, since these shift with distinctly different rates 
with temperature (typical values are 0.28 nm/K for the 
gain peak and 0.075 nm/K for the resonance, e.g. [20]). 
As considered in detail in [21] this mechanism controls 
the length scales of the transverse patterns emitted by 
the VCSEL. 

Figure [U shows the setup used for the experiments. The 
laser beam is collimated by a microscope objective with a 
numerical aperture of 0.8. VCSEL, cooling elements, and 
the objective are put into an air-tight box to avoid con- 
densation of water at low temperatures. The far-field of 
the laser emission is imaged onto a high-resolution 14-bit 
camera with a large charge-coupled device (CCD) chip. 
A half- wave plate and a linear polarizer are inserted into 
the beam path. The orientation of the polarizer defines 



the reference coordinate system by which the state of 
polarization is represented. Horizontal polarization is de- 
fined as 0°, angles are measured in counter-clockwise di- 
rection. The polarization is measured by taking far-field 
images for three settings of the polarization optics: hori- 
zontal {Ix), vertical (/^y), and diagonal orientation (145). 
For the circular component (Idrc) a quarter- wave plate 
(set to 45° with respect to the horizontal) is necessary. 
From this data, the spatial-resolved Stokes parameters 
are calculated: 
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where 5*0 represents the total intensity. Si the (normal- 
ized) amount of light polarized in x (positive ^i) re- 
spectively y (negative ^i) direction, S2 the (normalized) 
amount of light polarized along the diagonal direction 
(positive for 45°, negative for -45°) and S3 the (normal- 
ized) amount of circularly polarized light (the sign de- 
notes the direction of rotation). Using this set of Stokes 
parameters the degree of polarization (fractional polar- 
ization) p and the polarization direction cp can be calcu- 
lated: 



p = ^si 



SI 
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The fraction of circular polarization 5*3 was found to be 
of the order of 0.02, thus we assume linear polarization in 
the following. In this article we focus on the characteris- 
tics of two VCSELs, which illustrate the general behavior 
found in the experiments very well. 
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FIG. 2: (Color online) Near- and far- field images for de- 
vice #1 below and above threshold. Heat sink temperature 
T=0 °C and driving current 12.0 mA (a,d,g), 13.5 mA (b,e,h), 
and 15.6 mA (c,f,i). The color-code of the polarization di- 
rection distribution (lowest row) is shown by the bar on the 
lower right. These images correspond to the transition shown 
in Fig. [5l 



FIG. 3: (Color online) Near- and far- field images for device #2 
below and above threshold. Heat sink temperature T=0 °C 
and driving current 13.5 mA (a,d,g), 15.5 mA (b,e,h), and 
16.0 mA (c,f,i). The white spots in (a,b) result from debris 
on the neutral density filters. The color code is the same as 
in Fig. [2] These images correspond to the transition shown 
in Fig. [6] 



III. DISTRIBUTION OF SPONTANEOUS 
EMISSION AND THE TRANSITION THROUGH 
THRESHOLD 

We consider two nominally identical devices from the 
design and growth process, which show however rather 
different behavior. The near-field and far-field intensity 
distributions of the emission are shown in the figures [2] 
(device #1), resp. [3] (device #2). In both cases the heat 
sink temperature is T=0 °C, the threshold current for de- 
vice #1 is 15.2 mA, for device #2 15.6 mA. The first and 
second rows of each figure depict the intensity distribu- 
tions in grey-scale coding (black denoting the maximum 
intensity) of the near-field and far-field, respectively, the 
third shows the spatial frequency distribution of the po- 
larization direction in a cyclic color-code (red denotes 
±90°, green -45°, blue 0°, and yellow +45°, cf. to the 
color bar on the right of Fig. [2]). The columns show the 
emission below threshold (panels a, d, g; 12.0 mA for de- 
vice #1, 13.5 mA for device #2), slightly below threshold 
(panels b, e, h; 13.5 mA and 15.5 mA), and just above 
threshold (panels c, f, i; 15.6 mA and 16.0 mA). The 
optical axis is positioned in the center of each image. 

Below and above threshold, the emission has its max- 
imum at a well-defined wave number. This critical wave 
number is favored because it has the most favorable de- 



tuning properties as discussed above (see [2l[ for details 
about the dependence of the transverse wave numbers of 
the emission on the detuning). Even far below thresh- 
old the ring indicating the critical wave number is easily 
discernible. With current approaching threshold it nar- 
rows until at threshold the lasing modes develop from 
this ring. This is easily explained by the increase of the 
Finesse of the cavity if threshold is approached. 

Within this critical ring, far below threshold the max- 
ima of spontaneous emission is found at the diagonals in 
Fourier space, but moves close to the axes above thresh- 
old (i.e. with either small kx for device #1 or ky for de- 
vice #2). Just above threshold VCSEL #1 emits far-field 
patterns with two dominant Fourier components on the 
y-axis (Fig. [2](f)). The polarization of these components 
is in tendency orthogonal to their wave vector, which 
is shown in Fig. [^i), where the area of lasing emission 
is polarized horizontally (blue in the color-code). Be- 
low threshold the polarization direction is very different, 
parallel to the wave vector (see Fig.[2](g)). These general 
observations are also true for device #2. 

The validity of the '0°-rule' is illustrated and tested 
further in Fig. [H Here radial cuts through the polariza- 
tion distribution at 0°, 22.5°, 45°, 67.5°, and 90° with re- 
spect to the X-axis are shown for T=0 °C and 7=10 mA 
(i.e. far below threshold). Each curve starts at about 
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FIG. 4: (Color online) Radial cuts (at 0° (black solid line), 22.5° (red long-dashed line), 45° (green small- dashed line ), 67.5° 
(blue dotted line), and 90° (cyan dot-dashed line)) through the far- field distribution of the spontaneous emission for device 
#1. T=0 °C, 7=10 mA. (a) Polarization in dependence on the transverse wave number k± showing the validity of the 0°-rule 
for k± > 2/im~^. The thick lines indicate the polarization expected from the theory developed in Sect. II VI (b) Fractional 
polarization in dependence on k± for the same cuts as in (a). The part near k± ^ is shown in the inset. 



10° for k± ^ and switches asymptotically to the an- 
gle of the cut. For values of the transverse wave number 
above 1 /im~^ the polarization is clearly parallel to the 
wavenumber. The deviations amount up to about 10° 
at 1 /im~^ and decrease for higher wave numbers. The 
polarization state at /c^ ^ is interpreted to be selected 
by the cavity anisotropics. 

It is interesting to note that the curves of the polar- 
ization angle as well as the fractional polarization are 
continuous till the maximum value measured of 5 /im~^ 
though between 3.5 /im~^ and 5 /im~^ (depending on 
the direction of the wave vector) the cutoff condition for 
the transverse modes of the waveguide formed by the re- 
fractive index step (i.e. the side boundaries formed by 
oxidation) sets in. It is clearly visible in the center rows 
of Figs. [2] and [3] that the intensity is cut-off indeed. This 
indicates that the influence of the side boundaries on the 
polarization characteristics of below-threshold emission 
is very small. 

The fractional polarization, i.e. the amount of linearly 
polarized light, in dependence on the transverse wave 
number is shown in Fig.HJb), again for the same cuts. It 
increases monotonically with wavenumber. The graphs 
are more or less congruent which indicates the isotropic 
character of the phenomenon, showing its relative in- 
dependence from the principal axes of the intra-cavity 
anisotropy for large enough wavenumbers. The polariza- 
tion degree is small, but nonzero for k± = (see inset to 
Fig.HJb)), which also indicates the influence of the cavity 
anisotropics, because the DBRs are isotropic for k± = 0. 

In the following, we take a closer look at the changes in- 
volved in the transition from spontaneous to lasing emis- 
sion. This transition is illustrated in Fig. [5] for device 
#1 and Fig. [6] for device #2. Each figure shows in (a) 
the local (in Fourier space) intensity, in (b) the fractional 
polarization, and in (c) the local polarization orientation 
in dependence on the driving current. In addition the 



pump rate 

p = iLiM (3) 

is displayed on the upper x-axis of the diagrams. The 
inset in (a) shows the Fourier component for which these 
plots were made (indicated by the arrow). The black 
lines in the (a) panels show in both cases the typical 
behavior of a laser crossing threshold: A region of low 
emission intensity representing spontaneous emission and 
small slope is followed by a steeply increasing part indi- 
cating lasing emission. Threshold is extrapolated by a 
linear fit to this latter part and indicated by the dashed 
red line. The blue curves in Fig. [5] and [6] are calculated 
from equations (|9]), (p!2|) and will be discussed in the the- 
oretical section. 

The fractional polarization shown in (b) typically fol- 
lows the development of the intensity until it saturates 
at a maximum of 0.8 to 1.0. 

In part (c) of the figures the change of the polariza- 
tion direction with current is shown. Far below thresh- 
old the polarization angle is in good agreement with the 
0°-rule. With increasing current the scenario is differ- 
ent for the two lasers. For device #1, the polarization 
starts to change quite abruptly approximately 1.1 mA 
below threshold: In a current range of only 1 mA the an- 
gle changes to 0°, which corresponds to the "90-degree 
rule" . The polarization reaches the target state well be- 
low threshold (about 93% of Ith)- Note that the frac- 
tional polarization has a pronounced local minimum at 
the transition. 

Device #2 (Fig. El^c)) shows a rather gradual change 
of the polarization that starts already more than 5 mA 
(at 67 % of Ith) below threshold. The behavior is 
monotonous, there is no dip in the fractional polariza- 
tion. This behavior is typical, if the wavevector is not 
oriented along one of the two axis. Only in the latter 
case (Fig. [5]) an abrupt transition is found. The continu- 
ous transition is more typical for the devices under study. 
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FIG. 5: (Color online) Transition from spontaneous emission to lasing emission of the spot depicted in the inset in (a) for 
device #1. The black, solid curves are experimental data, the blue, dashed ones are calculated from Eqs. (|9))- (|12|) . (a) Local 
intensity in dependence on driving current. The threshold, indicated by the vertical red dashed line, is derived from a linear 
fit to the steep slope of the intensity, (b) Degree of polarization, (c) Local polarization direction. The change of polarization 
is evident about 1.5 mA below threshold. 



though, because the wavevector configuration depicted in 
Fig. [6^ is more typical [18|]. 

IV. THEORY AND DISCUSSION 
A. The Ginsburg-Landau equation 

In order to analyze the behavior described above 
we use a model for a broad-area VCSEL which ac- 
counts for its cavity structure including Bragg reflectors 
[13, For simplicity, we consider a spatially homo- 

geneous device with an infinite aperture. For this case 
the eigenmodes are plane transverse waves 'E\^^{x, t) = 
Yi{t) e-x.-p{i(kxX ^ kyy)} where E(x,?/,t) = {E^^Ey} is 
the slowly varying complex envelope of the field inside 
the cavity, = {kx^ky} is the transverse component of 
the wave vector. 

Many features of polarization selection at and slightly 
above threshold can be obtained by a linear stability 
analysis [l^, [2^. However, for the transition from be- 
low to above threshold it is of critical importance to take 
into account both spontaneous emission and the non- 
linear saturation. Therefore we will use here a nonlin- 
ear Ginsburg-Landau equation (GLE) with an additional 
term describing spontaneous emission (see Appendix El 
for the derivation). For the spatially homogeneous de- 
vice with infinite aperture the equation can be written 
for every transverse wave with complex amplitude E(t): 

E = -/^inE - A^outT(k^)E + z^](k^)E + FE- 

-/^oute(kx)/E + W. ^ ^ 

The field decay rate n = /^in + i^out results from the laser 
emission through the DBRs A^out (outcoupling losses) and 
intracavity losses A^in by scattering and absorption. The 
latter is isotropic (polarization independent), whereas the 
former is anisotropic. The anisotropy is described by the 
2x2-matrix T(k^), which represents polarization- and 



k^-dependent losses at the DBRs. In addition, the ma- 
trix T(kx) includes also the gain in the device (and hence 
has a component depending linearly on the driving cur- 
rent). l^(k^) represents diffraction in the cavity and in 
the DBRs, ^(k^) is a matrix describing the impact of 
the nonlinear saturation and / = E'^'E is the light inten- 
sity (t means the conjugate transpose). A more detailed 
description of T(kx), ^(kx) and ^(k^) is given in Ap- 
pendix [A) F is the intracavity anisotropy matrix, which 
in the basis of the main anisotropy axes is written as 
F = diag (7a + ^7p, —7a — ^Tp)^ where 7^ is the amplitude 
anisotropy (dichroism), 7^ is the phase anisotropy (bire- 
fringence) and diag ( • , • ) denotes a diagonal matrix with 
the corresponding entities on the diagonals. 

As indicated above, the other source of anisotropy is 
the reflection at the DBRs. The action of the DBRs 
can be described in terms of s- and p- waves [l^, HH, 
which are plane transverse waves with polarization cor- 
respondingly perpendicular and parallel to the direction 
of the transverse wave vector k^. In this basis, the ma- 
trix of reflection from the ith Bragg reflector is diago- 
nal Ri{\<it) = dia.g {Rsi{kt),Rpi{kt)}, i.e., pure s- and 
p- waves are reflected from the Bragg reflectors with- 
out mixing. The corresponding transmission matrices 
Ti{kt) = diag {Tsi(kt),Tpi{kt)}, z = 1,2 are also diag- 
onal in this basis. 

The spontaneous emission rate is described by the term 
W = W(kx,t) = ^(3spKj/Ti^{k^,t) in Eq. ^. Here 
^(k^,t) is a Langevin noise source (see Appendix [Al for 
details), j is the normalized current density, Ti is the 
population decay time, /3sp is the spontaneous emission 
factor (the fraction of spontaneous emission going into 
the given mode) , K is the Petermann excess quantum 
noise factor [26l], which takes into accounts a possible 
non-orthogonality of the modes leading to projection of 
the noise in other modes onto the selected one (26i] . 
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FIG. 6: (Color online) Transition from spontaneous emission to lasing emission of the spot depicted in the inset in (a) for 
device #2. (a) Local intensity in dependence on driving current, (b) Degree of polarization, (c) Local polarization direction. 
All denotations are analogous to Fig. O 



B. The coherence matrix 

In Eq. (|4j) the nonlinear term can be neglected for small 
current well below threshold, which results in the linear 
equation 

E = -/^inE - KontT{k^)E + in{k^)E + FE + W, (5) 
which can be solved directly: 

E = exp(6t) / exp(-6r)V^(r)dr, (6) 
Jo 

where O = —tZin — ^out^ + il] + F. 

If the field E is known, the Stokes parameters can be 
obtained from the coherence matrix J = (EE''') (here 
( • ) denotes an ensemble (and not time) averaging; note 
also the reverse order of multiplication in this definition 
compared to the definition of intensity, which results in 
a matrix instead of a scalar): 

Sj = tT{Ja,), (7) 

where aj are the Pauli matrices. In particular, the mean 
intensity (/) = can be obtained as (/) = tr J. 

By multiplying Eq. (|6]) from the left with its Hermite 
conjugate, performing averaging and then integration, we 
obtain 

J = -G-^Jw{l - exp(e^t)}, (8) 

where 0^ = + 0^ and Jw = (WW''') is the coher- 
ence matrix of spontaneous emission. In the derivation 
of Eq. ([8|) it is taken into account that the polarization 
components of W are delta-correlated and therefore Jw 
is proportional to the unit matrix Jwij = ^ijK jSgpj /2Ti 
(where J^j, i = 1,2 is a Kronecker delta), i.e., it com- 
mutes with all other matrices. It is also assumed that all 
functions of matrix arguments and 0^ commute. 

We will search for the statistically stationary solutions 
of Eq. (|8|) (i.e. those with a time-independent coherence 
matrix). For such solution being finite, the exponential 



term in Eq. ([8]) must decay. This is automatically ful- 
filled below threshold because is nothing but a linear 
stability matrix for the GLE without noise far below 
threshold near its non-lasing solution, and therefore all 
the eigenvalues of 0^ are less then zero. Hence, for the 
small current we obtain: 

J = -0-Vh^. (9) 

Eq. (O gives the coherence matrix inside the cavity. 
Because the Bragg reflectors transmit different polariza- 
tions differently, the coherence matrix Jo is different from 
J for an observer outside of the cavity: 

Jo = TJT\ (10) 

Let us turn our attention to the general case of the 
nonlinear Eq. (|4|). We can simplify the analysis by ne- 
glecting the joint fluctuation of the term /E and replace 
the intensity by its mean value (/) in this term. This 
approximation can be interpreted in the following way: 
We replace the original stochastic process described by 
Eq. (Hj) by a simpler, Gaussian one [27] with the same 
mean (E) = and, by its construction, with the same 
stationary mean intensity (/). We expect therefore that 
also the stationary coherence matrix J will not be signif- 
icantly altered by this approximation. Of course, the 
process described by the original GLE is not, strictly 
speaking, Gaussian, especially in the vicinity of threshold 
or polarization switchings [2^, [2^. However, as we will 
see later, this approximation fits rather good to the ex- 
perimental findings, allowing at the same time very con- 
structive analytical insight. By considering the resulting 
equation as a linear equation for the field, and proceeding 
as above we arrive at Eq. ([8]) with the modified matrix 

= -^in - ^outT -^iVt- noutQ{I) + r. (11) 

As before, we consider only the finite statistically station- 
ary solutions, for which the exponential term in Eq. ([8]) 
decays, and obtain again Eq. (|9|), with given now by 
Eq. (dH). 
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Taking the trace of Eq. ^ one obtains an implicit 
equation for the intensity (/): 



(/) =tr {(2Ki, 



-Kout(T + Tt)-(r + rt)+ 

+>^ont{i){g + g^)y' Jw}, 



final state instead of switching to it. In this case, the de- 
gree of polarization does not pass through zero and the 
polarization direction changes smoothly [see Fig. [6][b,c)]. 



(12) 



which is an equation of third order for (/). The intensity 
outside the cavity is then given by (Iq) = tr {TJT'''}. 
Eq. ([12]) has only one positive root, which is small 
((/) ~ 0) below threshold and grows asymptotically lin- 
early with current ((/) ~ j) above threshold. Below 
and at threshold it is however in rather good agreement 
with experiment (compare the black and blue curves in 
Fig. El^a) and Fig. [6Ka)). 



C. General features 

As stated before, the most significant difference be- 
tween the devices #1 and #2 is in the pattern formed 
above threshold. Whereas the below-threshold state 
obeys the "0-degree rule" for both devices (i.e. the po- 
larization direction is parallel to k^), the polarization 
of the above threshold state deviates significantly from 
the "90-degree rule" for device #2 (around 30°). In [3 
it was shown that for high enough k± the polarization 
direction is aligned mainly to the transverse side bound- 
aries of the device rather than to the anisotropics (either 
intra-cavity or Bragg-induced ones). The validity of the 
"90-degree rule" above threshold depends therefore on 
the position of the spots in the far field with respect to 
the directions of the side boundaries. For a device with 
the boundaries parallel to the x and y coordinate axes 
the "90-degree rule" is satisfied when the spots are lo- 
cated close to the x or y axes, whereas for spots away 
from the axes the polarization direction deviates from 
the "90-degree rule". Note that for the purpose of this 
paper it is not important why some devices emit in a 
specific wave vector configuration and some in another. 
We are only exploring the consequence of a given wave 
vector configuration on the polarization behavior. 

The difference between these two situations can be de- 
picted in the space of Stokes parameters (Fig. [7j). Be- 
cause the polarization is always linear, 5*3 = and we 
use only its two-dimensional subspace {Si^ 82)- One can 
see that - if the initial state E^^ and the final state E^/ 
are orthogonal to each other- 5*2 experiences a zero cross- 
ing during the transition (red cross in Fig.[7j). According 
to Eq. ([2|) this means that the degree of polarization is 
also zero at the crossing point. The polarization direc- 
tion retains only two values during this transition. Hence 
the transition, which appears at the point Si =82 = 0, 
is an abrupt switching between two discrete polarization 
states. Both the transition of the degree of polarization 
through zero and the abrupt switching can be seen in 
Fig. [5](b,c). On the other hand, when the final state Eg/ 
is not orthogonal to the initial one E^^, the polarization 
follows a path avoiding the origin, i.e., it rotates to the 




FIG. 7: (Color online) Different types of transitions from 
non-lasing to lasing state in the space of Stokes parameters 
{Si, 5*2) : Abrupt transition (blue straight dashed arrow from 
Eiai to Ea/, corresponds to device #1) and smooth tran- 
sition (green curved dashed arrow from Es^ to E^/, corre- 
sponds to device #2). The initial polarization state for the 
abrupt transition Eiai is orthogonally polarized to the final 
one Ea/, whereas for the initial Es^ and final E^/ states for 
the smooth transition this is not the case. For convenience, 
the polarizations directions for different quadrants of (5*1, 6*2)- 
plane are shown by red arrows in the black boxes. The red 
cross in the center corresponds to a fully unpolarized state 
{p = 0) whereas the black circle shows the fully polarized 
state (p = 1). 

We will refer to the first case as an "abrupt" transi- 
tion, and to the second case as a "smooth" one. One can 
see a remote analogy to the Ising and Bloch transitions 
[30], which also represent abrupt "switching" or smooth 
"rotating" behavior. The role of the vector of magne- 
tization in this case is played by the Stokes parameters. 
One should note however that Ising and Bloch transitions 
occur in space, i.e., between energetically equivalent spa- 
tially separated states (at constant external parameters), 
whereas here the transition takes place in dependence on 
an external parameter (current), i.e. in the parameter 
space. 

The above mentioned difference between the abrupt 
and smooth transitions can be expressed more directly in 
terms of the coherence matrix. For the abrupt transition 
the coherence matrix is diagonal in one and the same co- 
ordinate basis during the whole transition. Only two po- 
larization directions are possible, depending on which di- 
agonal element Jn or J22 is larger. The point represent- 
ing the polarization state in Fig. 8 can move only along a 
straight line passing through zero (^i = 0,5*2 = 0). At 
the point of transition Jn = J22, i.e. the light is unpolar- 
ized, and the polarization direction changes abruptly. In 
the case of a smooth transition the situation is different. 
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Of course, the coherence matrix is Hermitian and there- 
fore there is always a Cartesian coordinate basis in which 
it is diagonal. However, this basis is changing during the 
transition. That reflects the existence of several com- 
peting mechanisms (DBRs, side boundaries, intra-cavity 
anisotropics), each of them having its own principal axes. 
Their mutual influence is changing with a chance of pa- 
rameter (here current). In this case, the Stokes parame- 
ters do not vary along a straight line and can avoid the 
zero crossing. It should be noted also that the origin 
{Si = 0, ^2 =0) is a special point in the sense that the 
coherence matrix in this point is proportional to the iden- 
tity matrix and therefore is diagonal in every coordinate 
basis. 



D. Detailed discussion 

The polarization, intensity and fractional polarization 
obtained from Eq. (p!Q|) . Eq. ([9|) and Eq. ([12]) are shown 
in Fig. m O and [6] in comparison to experimental data. 
The parameters of the active layer and the cavity used 
for calculations are a = 3, 7 = 10^ ns~^, 7^ = 0.8 ns~^, 
7p = 0. The DBRs, consisting of 31 (top mirror) and 47 
(bottom mirror) A/4 layers of material with alternating 
refractive index ni = 3.46, 712 = 3.093 and the trans- 
parency current Jtr = 7 mA were assumed. The effective 
round-trip time r = 33 fs includes also the effects of the 
dispersion in the DBRs and in the cavity [21]. The out- 
coupling loses for these parameters are /^out = 17.4 ns~^. 
The best coincidence with the experimental results ap- 
pears for /^in = 65 ns~^, which is in rather good agree- 
ment with estimations based on losses in the p-DBR and 
typical gain values of GaAs quantum wells [3l|, [32| . 

In the framework of the theory presented here only an 
infinite device can be considered. Therefore, the final 
lasing state satisfies always the "90-degree rule" [I^ • As 
the laser approaches threshold, the influence of the side 
boundaries start to play an important role [18], and the 
polarization may not be perpendicular to anymore. 
For the device #1 this deviation is very small but it is 
not negligible anymore for the device #2. We take this 
into account in our model by introducing an artificial 
rotation of the main axes of the matrix for the device 
#2, so it is becomes diagonal not in the basis of s- and 
p-waves but in another, rotated one. On the other hand, 
we keep the transition matrix T the same (i. e. diagonal 
in the basis of s- and p-waves). This makes the resulting 
coherence matrix non-diagonal. The angle of rotation is 
25°, which corresponds approximately to the angle of the 
wavevectors visible in the inset of Fig. [3](a). 

Although the spontaneous emission factor is rather 
small {/3sp ~ 10~^) for VCSELs with a transverse size 
around 40 fim [33|, it can be significantly enhanced by 
the Peter mann excess factor. The value of K depends 
extremely strongly on the inhomogeneities and imper- 
fections in the construction of a particular device. The 
best results in comparison to experiment give the values 



Kf3sp = 6 X 10-"^ for the device #1 and Kf3sp = 10''^ for 
#2. In [18,] it is shown that the four spots at the corners 
of a rectangular which form the dominant Fourier peaks 
of the spatial structures (see inset of Fig. [6l^a)) can form 
the eigenmode of the transverse waveguide but are not 
simultaneously an eigenmode of the reflection operator 
of the DBR. Thus the reflection couples many transverse 
wavevectors . This may be the the origin of the rather 
large Petermann factor K, especially for device #2. In 
addition, non-orthogonality of the modes might originate 
from inhomogeneities of the structure and current distri- 
bution (which is clearly visible in the intensity distribu- 
tions in Fig. [2ja-c), Fig. [3]^a-c), for both device #1 and 
#2). 

With this choice of parameters, a very good agree- 
ment between experiment and theory is obtained for the 
development of the fractional polarization and the lo- 
cal polarization direction versus current for both devices 
(Figs. [5l and [6|) as well as for the dependence on the 
transverse wave vector (Fig. S]). The data reflect the de- 
gree of abruptness of the transition (Fig. [5t vs. Fig. [6t) 
and the monotonous vs. non- monotonous development 
of fractional polarization (Fig. [6)3 vs. Fig. [Sb). The in- 
crease of fractional polarization and the convergence to- 
wards the "90°— rule" with increasing wave number is 
due to the increasing anisotropy between s- and p-waves, 
of course (see the blue line in Fig. [8^1). 

The case of the abrupt transition allows more analyt- 
ical insight. Let us suppose for the sake of clarity that 
the isotropic intra-cavity anisotropy is diagonal in the 
representation of s- and p-waves. In this situation, all 
the matrices in Eq. (|9]) and Eq. (p!Q|) are diagonal in this 
representation. Therefore the coherence matrices J, Jo 
are also always diagonal. In this case Eq. ([9]) and Eq. (p!Q|) 
are decoupled into independent equations for the diago- 
nal elements: 

Kpspj 

2Ti{/^in + A>:outKe(Tii) + Re(rii) + /^out(^)Ke(fc)} ' 

Joii — Jii-) (1^) 

where i = 1 corresponds to the DBR s-wave whereas 
i = 2 to the p-wave. 

Well below threshold the denominator in Eq. (p!3|) is 
positive and approximately the same for both s- and p- 
waves (as it will be discussed below), and the light out- 
side of the cavity is slightly p-polarized (Joii < J 022) due 
to the filtering by the Bragg reflector. As the current 
approaches threshold the denominator for the s-wave in 
Eq. (p!3|) tends to zero whereas the one for the p-wave re- 
mains positive, which provides superiority for the s-wave 
strongly overcoming the opposite difference in transmis- 
sion. Obviously, at some point the transition between 
s- and p-polarization occurs. At the point of transi- 
tion J oil = Jo22 (and therefore Si = 0, ^2 = 0), i.e., 
the output is unpolarized (see Fig. [71 blue straight ar- 
row, and Fig. [Sfb)). Because the coherence matrix is 
diagonal, the polarization direction can have only two 
values, corresponding to either Jon < J022 ("0-degree 
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FIG. 8: (Color online) (a) The amplitude (blue, solid line) and phase ^hp (red, dashed line) anisotropy induced by the 
Bragg reflector, (b) The degree of polarization for the parameters of device #1 (see Fig. |4| and ky — ^ according to theoretical 
calculations. The polarization degree inside the cavity is indicated by the red dashed line, outside the cavity by the blue 
solid line. The corresponding experimental curve is depicted by the black line. The extra-cavity polarization in the "pure 
filtering limit" (i.e. assuming completely unpolarized light inside the cavity) is indicated by the green dot-dashed line, (c) The 
polarization direction inside (red dashed lines) and outside (blue solid lines) the cavity for the intra-cavity anisotropy favoring 
x-polarization (thick lines, as in panel b) or an anisotropy favoring a polarization direction of 10° (thin lines). 



rule") or Joii > Jo22 ("90-degree rule"). Therefore, at 
the transition point, the polarization direction changes 
abruptly and is constant in all other points above and 
below threshold (see Fig. [5jc)). 

Let us now consider the behavior of the extracav- 
ity polarization far below threshold (when the intensity- 
dependent term in Eq. ([13]) can be neglected). If we 
suppose that the intra-cavity losses K/iii are much larger 
than the ones through outcoupling /^out and than the in- 
tracavity anisotropy Tu we obtain: 

_ Kpspj t.._|^.|2t. 

In this case the output polarization is governed by the 
transmission through the DBR |T^^p, and the laser cav- 
ity works as a simple filter for almost unpolarized intra- 
cavity radiation. Because \Ts\^ > this results in the 
"0-degree rule" due to the dominance of p-polarization 
direction in transmission. In this case (we will call it 
"pure filtering case" here), the degree of polarization is 
determined only by the transmitting properties of the 
Bragg reflector (see Fig.[8fb), green line). The difference 
in transmission between s- and p-waves increases with 
k_\_ approximately quadratically, and the degree of polar- 
ization in this case reflects this dependence according to 
Eq. (dH). 

For the opposite case hZout ^ ^in and assuming a diag- 
onal coherence matrix, we obtain 

Jou-\Tu\VRe{Tu), (15) 

instead of Eq. Because for small current T^^ ~ 

1 — \Rii\'^ = |Tiip, the light outside the cavity is com- 
pletely unpolarized in this case. This is in agreement with 
the energy conservation principle, since the numbers of 
intracavity photons in both polarizations are increased 
by the spontaneous emission equally and the polariza- 
tions do not mix. Therefore, in the stationary case in 
the absence of intracavity losses the energy escaping the 
cavity must be also equal for both polarizations. 



In general, the extracavity polarization degree, defined 
by Eq. ([13]), lies in between these two limiting cases 
(Fig. [8](b), thick, solid blue line). As the Bragg reflec- 
tion is isotropic for k± = {Ta I^^^q ~ -'- ^ ~ l?^), 
the polarization for is defined only by the intra- 

cavity amplitude anisotropy (and does not depend on 
7p). The best agreement with the experimental value 
of the degree of polarization for = is achieved for 
7a = 0.8 ns~^, a reasonable number in line with typical 
observations in small-area VCSELs [H, [35[. With in- 
creasing transverse wavenumber k±, the relative impor- 
tance of the anisotropy induced by T(kx) increases. The 
amplitude anisotropy Jba0^±) = ^out(R'e(Tii(k^)) — 
Re(T22(kx)))/2 and the phase anisotropy 76p(k^) = 
Kout(Im(Tn(ki)) -Im(T22(ki)))/2, induced by T(ki) 
are shown in Fig. (Hj^a). In analogy to the intra-cavity 
anisotropy, only 75^ plays a role in the determining of the 
coherence matrix. As one can see, 75^ increases approxi- 
mately quadratically with k± . The degree of polarization 
inside the cavity (see Fig.jSj^b), thick, dashed red line) is 
influenced more strongly by the intracavity anisotropy for 
small k± (below ~ 1.5 /im~^) and by the Bragg-induced 
anisotropy for larger k±. But it remains relatively small 
for all k± , and the polarization degree outside of the cav- 
ity is therfore quite strongly determined by the Bragg 
filtering mechanism (see Fig.jSj^b), thick, solid blue line). 

Now, let us consider the behavior of the polarization of 
intra-cavity field in dependence on k^ far below thresh- 
old in more detail. As an example, a cut displaying 
the polarization angle along the kx axis, assuming the 
anisotropy also being directed along the x-axis is shown 
in Fig. Fig. [8](c). In this case, the intracavity and ex- 
tracavity light for k± ^ is x-polarized. However, for 
high enough k± (above ^ 1.5 /im~^) the DBR-induced 
anisotropy overcomes the intracavity one and the intra- 
cavity light becomes weakly polarized in the direction 
perpendicular to k^ (see Fig. |8]^c), thick red, dashed 
line). This degree of polarization is small and canceled 
when the light is transmitted through the Bragg reflector. 
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Hence the polarization of the hght outside of the cavity 
is stih determined by the transmission, i.e., is parallel to 
kx, as it mentioned above (see Fig. [81(c), blue lines). 

Hence, we also encounter a polarization transition in 
the intracavity field, if we consider k± as a parameter 
instead of the current. This transition is abrupt (see 
Fig. (He), thick red, dashed line) and the fractional po- 
larization vanishes at the transition point (see Fig. [Hl^b), 
thin red, dashed line). It can be understood in a way 
fully analogous to the abrupt transition appearing with 
the change of current. It should be noted, that this tran- 
sition is not apparent outside the cavity (thick blue line 
in Fig. He)). 

This behavior becomes slightly more complex, if one 
analyzes the polarization along a cut in Fourier space 
whose direction does not coincide with the preferred 
intra-cavity anisotropy axis. In that case, the change of 
polarization angle vs. wavenumber is smooth (Fig. [8l^c), 
thin red line) and also detectable outside the cavity 
(Fig. [81(c), thin blue line). This is backed up by the exper- 
imental curves displayed in Fig. H^, where the transition 
is most notable for the cyan curve, which corresponds 
to a cut along 90°; i.e., the competition between intra- 
cavity anisotropy and Bragg anisotropy is more apparent 
in the extra-cavity polarization, if the deviation between 
wave vector and anisotropy angle increases. 

V. CONCLUSION 

For highly divergent emission in wide aperture VC- 
SELs the polarization direction is determined by the 
properties of the distributed Bragg reflectors (DBRs), 
and (close to threshold) by the transverse structure of 
the cavity. Far below threshold the polarization is in- 
dependent from the cavity structure and the p-wave of 
the DBRs (which has the polarization direction parallel 
to k^) prevails in the emission because the almost un- 
polarized intra-cavity light is filtered by the transmission 
through the DBR, which is higher for p-waves than for 
s- waves. In the simplest case, when the above thresh- 
old polarization direction is not strongly influenced by 
the transverse boundaries, an abrupt switch of polar- 
ization direction occurs as the current increases towards 
the laser threshold: the s-wave inside the cavity start to 
prevail because the reflectivity and therefore the quality 
factor of the cavity is higher compared to the one for p- 
wave. At the point of transition the filtering effect of the 
transmission through the DBR can not compensate the 
intra-cavity polarization anymore, and the polarization 
changes its direction from the one parallel to k^ to the 
perpendicular one. This abrupt change of polarization 
direction is accomplished by the passing through zero of 



the degree of polarization of the light outside of the cav- 
ity. 

On the other hand, if the above-threshold polarization 
does not coincide with the DBR s-mode, the transition 
is qualitatively different. In this case, the polarization 
changes smoothly and the degree of polarization is signifi- 
cantly nonzero during all the transition. As was shown in 
[18], the deviation of the above-threshold polarization di- 
rection from the one dictated by the s-wave of the DBRs 
can be due to the influence of the side boundaries of the 
cavity. 

The difference between two types of transition can be 
explained in the terms of coherence matrix. The abrupt 
transition occurs when the coherence matrix of the light 
outside the cavity is diagonal (in one and the same ba- 
sis) during the whole transition. In this case, only two 
polarization directions are possible. If the different mech- 
anisms influencing the polarization have different direc- 
tions, the resulting coherence matrix is non-diagonal and 
arbitrary polarization direction is possible. 

We note an analogy between the abrupt and smooth 
transitions described here to Ising and Bloch transitions 
between two equivalent states in ferromagnetics with re- 
spect to a "switching" vs. "rotation" behavior. However, 
the polarization transition in this article occurs in param- 
eter space (between energetically inequivalent states) in 
contrast to the usual Bloch and Ising ones. 
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APPENDIX A: THE GINSBURG-LANDAU 
EQUATION 

1. The initial equations 

Polarization phenomena in VCSEL are often modeled 
using an equations for the intracavity field E and the 
total carrier density D as well as the population difference 
d between sub-bands with opposite carrier spin B, . 
We start from the nonlinear equations (18) of [13] for 
the normalized complex- valued envelope E(t, k^), of the 
optical field for given k^ and carrier population variables 
D, d: 
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E = -A^outME-inE-moutQ^E + A€out(l + ice)G(^E)+W, (Al) 
b = -7i{L>- j + Im[(i-a)E*/:(AE)]}, (A2) 
d = -7,(i-7iRe[(i-Q^)E*£(AE0]. (A3) 



Here a is the hne width enhancement factor, j is the 
normahzed current density, tZout is the outcouphng losses 
(see below), 7i = 1/Ti and 7^ are the decay rates of D 
and d, E' = {Ey,-E^} and A = (4^*^). The lin- 
ear operators M(k_L), 17(k_L), G(k_L), jC(k±) describe the 
losses, diffraction and gain in the cavity, and are 2x2 
kx-dependent matrices acting on the vector field E(kx). 

In the present paper we assume that r^(k^) = 
—k'^v/kol^ + (<si(k_L) + S2(k_L)) /r is a matrix describ- 
ing the dispersion relation given by the cavity reso- 
nance condition for s- and p- waves of DBR. Here v 
is the speed of light in the cavity, ko is the longitudi- 
nal part of the wavevector, r is the cavity round-trip 
time, \^ = diag (1,1) (diag ( • , • ) is here a diagonal ma- 
trix with corresponding entities on the diagonals), and 
si, S2 are the matrices, describing diffraction of the light 
in the DBRs. In s-p-representation they can be written 
as Si = diag (sgi, 5pi), where i = 1,2 and Ssi^Spi are the 
phase shift for s- and p- waves. 

Using these assumptions, M(k_L) and G(k_L) 
can be written in terms of propagation ma- 
trices Fi(k_L) = (pf)V2i^i(k_L)exp(-isi(k_L)), 
F2(k^) = (pf)V2i^2(k^)exp(-i52(k^)), F = 
F1F2: M(k^) = (l-F(k^))/Mo,G(k^) = 
(1 + Fi(k^) + F2(kx) + F(k^)) Cn/Go. M and 

G are normalized by the constants Mq and Go 
in such a way, that Mii(O) = 1, Gii(O) = 1. 
Here i^(kt) is an operator describing the refiection 
from the Bragg mirrors, represented by matrices 
Rm = Rmij{m l,2,i = 1,2, j 2,2). Kout = 

(1 - |i?lll(0)||i?21l(0)|)/T = (1 - |i?122(0)||i?222(0)|)/r 

is the outcouphng losses for zero transverse mode, p 
describes all intracavity losses in the whole system 
(which are not due to outcouphng) and V is the intra- 
cavity anisotropy matrix, which in the basis of the main 
anisotropy axes is written as 

p _ Aexp(7.r+Z7pr) « ) (A4) 

where 7^ is the amplitude anisotropy and 7^ is the bire- 
fringence. The infiuence of the gain contour line C{}<i±_) is 
given by the expression C=l/ {l ^ {{3 — Q) / ^)'^) where 
L is the cavity length and 7 is the material polarization 
decay rate. 

The spontaneous emission is described by the term 
W = KPspJiD^ in the approximation of zero inter- 
sub-band population difference. Here, K is the Peter- 
mann excess quantum noise factor and /3sp is the spon- 
taneous emission factor. ^ is the Langevin noise source 
with zero mean and correlation in (x, 7/) -space and circu- 



lar wave basis ^{x,y,t) = {^^{x,y,t), ^-{x,y,t)}: 

(e± {x, y, t)e± (x^ y\ t')) = 23{t - t')3{x - x')3{y - y'), 
{^±{x,y,t)^^{x',y',t')) = 0. (A5) 

Performing the transverse Fourier transform and trans- 
forming into a basis of linear (orthonormal) polariza- 
tion with arbitrary directions of the axes 1 and 2 
one obtains the analogous equation for ^{kx^ky^t) = 
{^1 {kx , , t) , ^2 {kx ,ky,t)}: 

{i^{kxMUK.K^t')) = 23ij3{t-t')x 

x3{kx-k',)3{ky-k'y). ^ ^ 

Therefore, the noise is also correlated in the k^- 
representation and in an arbitrary orthogonal polariza- 
tion basis. The noise terms in the equations for D and d 
are neglected in the present consideration. 



2. The derivation of the Ginsburg-Landau equation 
(GLE) for the field 

Here we obtain the lowest-order nonlinear equations 
for the field resulting from the above mentioned nonlin- 
ear equations. We take into account that near lasing 
threshold the resulting linear operator acting on E, D, d 
has a block-diagonal form with the part acting on E not 
being coupled to the carrier part. The eigenvalues stem- 
ming from the carrier related part are always strongly 
negative. In addition, the solution with (i = is always 
possible. Then, it is possible to adiabatically eliminate 
D and obtain a complex equation for E only. 

The solution of I) = is given then hy D = j / (1 -h/>C), 
where / = |Ep. For small intensity / one can write it as 
D ~ j(l — IC). Substituting this into the equation for 
the field, we obtain: 

E = -A^outME - i^E - iA^outceE -h /^out(l + ia)jGCE- 

-/^out(l + ic^)jG/:2/E + W. 

(A7) 

Because the anisotropy and the intracavity losses are 
small compared to 1 , the corresponding terms can be de- 
composed as p = exp(— /^inr) ^ 1 — /^inT, F ^ 1 + rr with 
F = diag (7a + i7p, -7a - ilp)- Considering /^inT, /^out^ 
and Fr as small parameters and neglecting these terms 
starting from the first order in G and from the second or- 
der in M, and introducing the matrices Q = {l+ia)jGC^ ^ 
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G 



-- M 
G 



(1 + ia)jGC (where M = M 



=i,r=i ' 



=i,r: 



^-^ ) we obtain the resulting GLE (j3|). 



tion of a smah intensity can be written as W(kx,t) = 
^y KPsplij^O^±^t). Here we neglected the second order 
term in decomposition of stationary value of D into se- 



The spontaneous emission term in the approxima- ries. 
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